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SUPREMUM DISTRIBUTION OF BESSEL PROCESS OF DRIFTING 

BROWNIAN MOTION 


A. PYC, G. SERAFIN AND T. ZAK 

Abstract. Let + /it) be a three-dimensional Brownian motion with drift /i, 

starting at the origin. Then Xt = \\{B^\B^\Bf'^ + its distance from the starting point, 
is a diffusion with many applications. We investigate the distribution of the supremum of (Xt), 
give an infinite-series formula for its density and an exact estimate by elementary functions. 


1. INTRODUCTION 


In his famous paper ([9], see also |8], page 436) David Williams showed how one can decom¬ 
pose the paths of a transient one-dimensional diffusion at its maximum (or minimum). One of 
the best-known examples of such decomposition is that of B{t) -|- fit, Brownian motion with a 
positive drift p, as a Brownian motion with a negative drift B{t) — fit and a diffusion with a 
generator 


( 1 . 1 ) 


A. 


1 d? 

2 dx^ 


-|- fi coth(/ia;) 


d 

dx 


Namely, for /i > 0 let B{t) —fit be a Brownian motion with constant drift —fi, Zt a diffusion with 
generator (HH), started at zero, and let 7 be a random variable with exponential distribution 
with parameter Assume that Bt, Zt and 7, defined on the same probability space, are 
independent. Define r = inf{f > 0 : Bt — fit = —7}. Then the following process 


( Bt — fit, 0 <t < T, 

\ Zt-r - 7 , t>T, 


is Bt + fit, a Brownian motion with positive drift fi. Williams also showed that Zt can be viewed 
as a Brownian motion with drift conditioned to hit infinity before zero. 


Later on a process with generator ma, known as a hyperbolic Bessel proces (0, p. 357) or 
as a Bessel process of drifting Brownian motion and denoted BES(3,/i) (|6], |3]), appeared also 
in many papers as a drifting Brownian motion conditioned to not hit zero. Recently, it was 
proved in pQ that it can be obtained as a deterministic involution of Brownian motion with drift 
fi. li fi = 1 then (Zt) is a radial part of a hyperbolic Brownian motion in three-dimensional 
hyperbolic space. 


The transition density function of {Zt) is well-known, compare [6] or |3] but to our best 
knowledge, distribiution of different functionals of this process have not been investigated yet. 
In this paper we investigate process {Zt) killed on exiting interval (0, tq) and give a formula 
describing distribution of Mt = supg<jZs, the supremum of the process {Zt). Because the 
formula is given as an infinite series, we give its exact estimate using elementary functions. 
Moreover, our method of estimation applied to a function ssy{v,t) used in a handbook by 
Borodin and Salminen give very precise estimate of this function (cf. Remark after Theorem 
9). 
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2. Transition density of a Bessel process of drifting Brownian motion, 

KILLED ON EXITING INTERVAL (0, Tq) 


Let Z{t) be a diffusion on (0, cxd) generated by operator (1.1) with p > 0. The speed measure 
of this diffusion is equal to m{dy) = sinh^(/i?/) dy, hence with respect to m{dy), the transition 
density of Z{t) has the following form (cf. [B]): 


( 2 . 2 ) 


p{t;x,y) = 


l{2t) _ g-(3/+a:)2/(2t) j 

\/27rt sinh(/ia;) sinh(/i?/) 


By (Z[°) we will denote the process killed on exiting (0, ro). For positive y, the process starting 
from X > 0 cannot reach zero and drifts to inhnity so that almost all trajectories will be killed 
at tq. Even if (Z^) starts from zero, with probability one it will never visit zero again. 

Transition density p‘^°{t-,x,y) of (Z[°), with respect to m{dy), is a solution of the following 
Dirichlet problem: 


f t > 0, X e (0,ro), ye (0,ro), 

(2.3) < p^o(t;x,ro) = 0 t > 0, x G (0, Tq), 

p^°{t;x,y) = S^{y) x e (0,ro), y G (0,ro). 


Because of killing, p^°{t;x,ro) = 0 for all t > 0. Moreover, if /r > 0 then by 02 .21) we 
have limsup^^op(t; X, y) < oo for all t, x > 0 and, because p^°{t;x,y) < p{t;x,y), hence 
\ims\ipy^QP^°{t;x,y) < oo. 

We will use the separation variable technique used in mathematical physics. Suppose that 


p^^{t-x,y) = Y{y)T{t). 

Then, by 01.ip . the hrst equation of the system 02.31) takes on the form 

(2.4) Y{y)r{t) = fly"(a) + acoth(ra)y'(a)j T(t). 

We add |(A^ + pi?)T{t)Y{y) to both sides of 02.4|) and get 

Y ( y )\ T '( t ) + 1(a2 + P ) T ( t )] = r(i)[ly"(») + ycoHm ) Y '( y ) + i(A^ + a")y(9)]. 

The solution of the differential equation T'{t) + A(A^ + pi?)T{t) = 0 is T{t) = cic \ If in 
the second equation 

Y''{y) + 2/icoth(/i?/)y'(?/) + (A^ + p^)Y {y) = 0 


we substitute Y {y) 


u{y)/ sinh(/i 2 /), we get 


1 

sinh py 


{u"{y) + XMv)) 


0 . 


By the above discussion we know that limsup^^^g ^(l/) < '^iv) = sinh(/i|/), hence 

\miy^Qu{y) = 0. This implies that with the boundary conditions m(0) = uiyo) = 0 equation 
u''{y) + X^uiy) = 0 is a special case of the general Sturm - Liouville problem [5]. It is well 
known (cf. [5], Theorem 4.1 and Exercise 4.2, page 337) that this Sturm - Liouville problem for 
u{y) has a solution if and only ifA = ^,n=l,2 ,... and this solution (up to a multiplicative 
constant) is given by Un{y) = sin(?7,7r|//ro). Thus we may expand p^°(t; x, y) as 


n=l 


an{x) , /nvri/AA 

sinh{py f"^\ To JJ 


exp 


/ -(nV/rS+ 


(2.5) 
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sin{^) 

. ^ rn ' 


To determine coefficient a„(a;) let us multiply the equation fl2.5p by Qfcand integrate 
the product over (0,ro) with respect to the measure sinh^(/ir/)(ir/: 


( 2 . 6 ) 


Io° y)ak{x) sinh^ {yy)dy 


_ rro 

Jo 2_-/n=l 


sin f sin f 


sinh^(/iy) 


Sin 


\ ro J \ ro J 


e^p ^ 


smlr^ {ny)dy. 


Observe that 


o sin(^) sin(^) m 

an{x)ak{x) sinh^(p,y) . ,x . , -. dr = a^{x)—6'^, 


sinh(/ii/) sinh(/ir/) 


where 5(1 is the Kronecker delta. Hence the right-hand side of 02.61) is equal to 


^0 2 / \ 

—a„(a;)exp 




2r§ 


Now, let t —?• 0. We have exp using the third condition from 02.3p for 

the left-hand side of 02.6p . we get 


afc(x) sin 


( knx \ 

ro J 


sinh(/ra;) 


^0 2 / \ 

= 


hence 


ak{x) = 


2 sin 


/ knx 

V ro 


ro sinh(/ix) 

To sum up, we have just proved the following theorem. 

Theorem 1. Transition density (with respect to the measure sinh^(/ir/) dy) of the Bessel process 
of drifting Brownian motion, starting from the point x G (0,ro) and killed at ro, is given by the 
following formula 


(2.7) 


fit, x,y)^Yi 


n=l 


2 sin ( ^ 

ro 


sm 


/ TlTTy 

\~ J 


ro sinh(/ia;) sinh(/ir/) 


exp 


[rrirlr\ -|- pB)t 


We will give another representation of the transition density, using the Poisson summation 
formula (cf. 13.4 in [ID]): for any function g absolutely integrable on (— cx), cxo) 


(2.8) 


CXD CXD 1*00 

9i.n) =Y d{x) 

7- J —OO 


—27Tikx 


dx. 


k=—oo ‘ 


First, observe that series (ETP can be written as 


p 


TO 


if]x,y) = 7; Y 


2 sin 


(^) (1?) 

ro sinh(/ix) sinh(/i?/) 


exp 


lr\ -I- /i^)f 


For the Poisson summation formula we need to compute the Fourier transform of any term 
of the above series, taking u as a variable. Put 2 ; instead of n, then we have to evaluate the 
following integral: 


sm 


Z'KX 


FO 


sm 


ZTxy 

Fo 
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In order to compute it, we use trigonometric identity e = cos(27rfc^) — ism{27rkz) and 
observe that the integral with sin(27rfcz) is zero, hence we have to compute only 

,in sin co.{2„kz)e=itdz. 



Now we use identity 2 sin ax sin bx = cos((a — b)x) — cos((a + b)x) and a formula (3. 
[1] to get: 


.2) from 


sin 


' —OO 


V ^0 


ro 


dz = 


Tq ( + -{y+x- 2 krf))^ 

e 2t — e 2t 


+ e 


V 

The Poisson summation formula fl2.8p gives us the following: 


-(y-x-2fcrn)^ -{y+x+2krQ)‘^ 

2t — g 2t 


Theorem 2. Transition density (with respect to the measure sinh^(/ij/) dy) of the Bessel process 
of drifting Brownian motion, starting from the point x G (0, ro) and killed at r^, is given by the 
following formula 


x,y) 


e 2 


\/27rf sinh(/ia;) sinh(/4|/) 


E 


— (y-x + 2krf)p —(y+x+2krQp 

e 2 t — e 2 t 


3. Exit time 


Let us consider Mt = supg<j Zg, the supremum of the Bessel process of drifting Brownian 
motion. The distribution of {Mt) is closely related to the distribution of the time when the 
process {Zt) exits the interval (0, ro). Recall that for /r > 0 process {Zt) exits (0, ro) at the 
point ro. For ro > 0 let us dehne = inf{s : Zg > ro}. Distribution of {Mt) and the survival 
probability of the killed process (E}°) are related by the following formula: 


¥^{Mt < ro) = P"(r,g > t) 


-ro 


p''^{t]x,y)sin\i^{yy)dy. 


Jo 

For hxed fv and ro series fl2.7p is uniformly convergent for x, y E [0, ro], because there 
holds inequality I I — and ^ exp(—n^7r^/(2rQ)) < oo. Thus we may 

integrate the series term by term. Integrating n-th term, we get (compare |1], formula 2.671) 




smh.{fiy)dy 


(_ 1)«+%r on sinh (yur 0 ) 

(n^TT^ + /i^rg) ’ 


so that we may write the following: 


Theorem 3. For t, r^ > 0 and x G (0, rg) the following formula holds 


F^{Mt < ro) = P"(r,g > t) = 


(3.9) 


E 

n=l 


(_l)n+i 2 ^^ sinh(/iro) sm{mrx/rg) 


smh{fix){n‘^7i‘^ + pBr, 


exp 


(nV/rg + ia^)t 


If we differentiate the above series term by term with respect to t we get 


71 sinh(yuro) 
sinh(/xa;)r| 


2 exp 


pft 

~2 


OO 


y {—l)^nsm{mrx/rg) exp 


n=l 


n^TT^ 
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Fixing e > 0 we can observe that the series of the derivatives is convergent nniformly for 
t G [e, oo) so that differentiation term by term is justified. Since G dt) = —^P*(rro > f), 

the exit time density is given by the following formula. 


Theorem 4. For fixed tq > 0, 0 < x < tq and any f > 0 

(3.10) P^(r„ G dt) = exp ( V(-l)’"+^nsin(n7rx/ro) exp 

smh(/ia;)r^ V ^ J 


n‘^71^ 


In a similar way as we did it for the density of the killed process, using the Poisson summation 
formula we can obtain another representation of the exit time density. Note that the series in 
formula fld.lOp can be written down in the following form: 


/ 2 2 \ 1 °° / „2 2 \ 
l)’^+^nsin(n7ra;/ro) exp ( — ^ tj = - ^ (—1 V sin(n7ra:/ro) exp ( —' 

In order to use 02.81) . we will compute the Fourier transform of its n-th term, taking z in place 
of n. First we take only cosine of exp{—27iikz), next we integrate by parts and hnally we use 
(0, formula 3.896.4) to get 


(3,11) 


2;sm 


ZTiy 

ro 


exp 


-z^ttH 


2rl 


exp{— 27 iikz)dz = 


\/2(7rf)3/2 _ 

Observe that 


{y + 2fcro) exp 


■{y + 2kro) 


sm 


2t 

mi{ro — x) 
ro 


+ (y — 2kro) exp 


■(y - 2kro) 
2t 


= (—l)"’’*'^ sin(? 7 , 7 rx/ro). 


hence putting y = tq — x m. 03.lip and using Poisson formula, we get the second representation 
of the exit time density. 

Theorem 5. For fixed fq > 0, 0 < x < fq and any t > 0 

TT^r/ sinh(uro)e“'^^*/^ ^ , , / —(ro — x + 2fcro)^ 

G dt) = . ^ Y1 i^o-x + 2kro) exp ' ^ 


sm 


h(/ix)\/^f3/2 


k=—oo 

4. Mean exit time 


2t 


Now we want to compute E*(rrQ) — the mean exit time of (Zt) from (0, ro). We will use the 
formula 

poo 

E"(r,J= / tF%Tr,edt), 


so that we need to compute the integral 


sinh(/iro) 


fOO —y?tl2 °° 


(ro — X + 2fcro) exp 


-(ro - X + 2/cro)' 


\/^sinh(px) Jo Vi 
Integrating a single term we use formula (3.471(15)) from |1] and get 

POO 2 


2t 


dt. 




-(ro + 2/cro - x)2 Vt\ V^ ^ i , o/ l^.7^ 

- — - YJ ^ exp(-/i|ro + 2kro - x\)dt, 
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which gives 


(4.12) IE"'(rr-o) = ^sinUjj^x) ^ [{ro - x + 2kro) exp{-fi\ro - x + 


2krn\)] = 


k=—oo 


sinh(/iro) 

/isinh(/ix) 


exp(/i(-ro + x - 2kro))] + 

OO 

[(ro - X - 2kro) exp(/i(ro - x - 2kro))] + (ro - x) exp(-/i(ro - x)) 

.k=l 


But 


and 


(ro - x)(e^("°-") + ^exp(-2/ifcro) = 


fc=i 


2ro(e^("-^°) -e'^(^°-"))f;fcexp(-2pfcro) = 


A:=l 


If we put them into fl4.12p . after some algebraic manipulation the formula for W{tb) can be 
simplified a lot. Namely, we get the following: 

Theorem 6. For any fixed tq > 0 and any starting point x G (0, ro) 


~ (’"o coth(pro) — xcoth(p,x)) 


5. Estimates 


Except for the last one, all the above formulas are given as series so that they are not 
convenient for computations or applications. In this section we give exact approximations by 
elementary functions of the transition density of the killed process, of the killing time and of 
the density of the distribution of the snpremum of the process. Notation f ^ g means that 
there exist two absolute constants Ci and C 2 snch that for all possible values of variables and 
parameters there holds Cif < g < C 2 /. 

For simplicity let us denote 7 '’°(t;x) = P* (r^p G dt). Recall that by Theorems 4 and 5 

7’'«(t,n = 


(5.13) 

(5.14) 




_ smh.{firQ)e ^ { 


X + 2nro) exp 


(rp—x+2nro)^ 
2i 


and, by Theorem 2, 


(5.15) 


g —M t /^^00 
sinh(/i3:) sinh(//y)\/27rt 2-^n=—oo 


exp 


(y—xfi-2nrQy 

2t 


exp 


{y+x+2nrQy 

2t 


We will start with estimates in the particular case ro = 1. We also make separate calculations 
for t G (0, i] and for t G [|, 00 ). 
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Theorem 7. Let tq = 1. For 0 < t ^ and 0 < a; < 1 


0,25 < 




1 x{l-x) Sinh/J -;,2^/2-(1-.r)2/r2U 
■^/^t3/2 t+x sinh(/xa:) 


< 4,02. 


Proof. First we consider the case 0 < a; < | and nse fl5.14p . It is enongh to estimate only the 
part of 7 ^(t,x) consisting of a series, mnltiplied by jg following qnantity 


M-x)^/{2t) °° 

/ =- ^ (2fc + 1 - a;)e-(2fc+i-^)V(2i). 

k=—oo 


X 


Observe that we can gronp terms of the series: k = 0 with k = —1, k = 1 with k = —2 and so 
on. In this way we get 


(5,16) l = T,<= 

k=0 

Now 


-2k(k+i-x)/t {2k + l-x)-{2k + l + x)e 2(2fc+iWt 


X 


a+.)/ = a+.) V ^ 

X 

k=0 

^ g- 2 fc(fc+i-x)A f2k + l-x-t + -{2k + 1)(1 - - (2A: + 1 + X + t)e-2(2fc+iwA . 

k=0 \ X ) 

Below we will estimate separately the term for /c = 0 and terms for k>l. The term for fc = 0 is 
eqnal to g{x,t) = 1 — x — t + |(1 — — (1+ a; + t)e“2^A. We will show that for 0 < x < 1/2 

and 0 < f < 1/4 there holds g{x,t) < 2. Indeed, a derivative 


dg{x,t) e 2^A 


(2x2(1 + x) - tx(x(e2"/* - 1) - 2) - t2(e2-A - i)) 


dx tx2 

is negative for x > 0 and f > 0, becanse _ x > ^ _|_ ^mX this inequality implies 
that in the formula for the derivative, the quantity in the bracket is negative. Hence g{x,t) < 
linia-^o g{x, t) = 2 — 2t and this is less than 2 for t > 0. 

In order to estimate terms for /c > 1, we use assumptions 0<t<|,0<x<| and inequality 
1 — < X, valid for x > 0 , and get 


{t + x)I <2 + (2k + l + ^{2k + l){l 


k=l 


^-2{2k+l)x 


/*)) < 


14 


2 + ^ (2A: + 1) (4A: + 3) < 2 + 21e"® + ^ < 


k=l 


k=2 


14 


2 + 21e-® + — V ne-^ < 2, 01. 


n=32 


To get an estimate from below, we will use inequality 1 — e ^ valid for x > —1. Using 

formula fIS.lOp . we get for 0 < t < ^ and 0 < x < |: 


fc =0 


—2k{k+l—x)/t 


-2x + {2k + 1 + x){l - e-2(2fc+iAA) 


> 


E' 

k=0 


,-2fc(fc+i-x)/f A A 1 

t + 2{2k + l)x 


X 

CXD 

k=0 


-2fc(fc+i-x)/t 2fc + 1 - ^ 


t + x 
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oo 

— Y1 + -) > 


t + X 


4^ 4(t + a;)’ 


so that {t + x) I > Both above estimates imply the following 

, I oo 

0, 25 < (2k+l- a;)e-(2^+i-^)V(2t) < 2 , Ql. 


If, for 0 < a; < |, we divide the middle term of the above inequality by (1 — x) G [|, 1), we 
must multiply its left and right-hand sides by, respectively, 1 and 2. This proves the theorem 
in this case. 

Now the proof for the case | < a; < 1 and 0 < t < |. Similarly as in the proof of the hrst 
case, we will examine the following quantity: 

(l—x)^/( 2 t) 00 ..00 

^ {2k+l- _ --- ^ (2A; + 1 - = 


-.00 OO 

_ 2 A; _ g-( 2 fc^- 2 fc(l- 3 :))/t'\ ^-{ 2 k^+ 2 k{l-x))/t ^ ^ ^ 

1 CC \ / 1 cc 

k=l k=—oo 

First we estimate the series denoted by 5; for 0 < 1 — a; < ^ and 0 < t < ^ 

OO OO 

B — 1 = 'Y^ _|_ g-(2fc^-2fc(i-x))/t j ^ 'Y^ -f < 


(e-«“ + e-““ 


S3 + 4+E2e-‘<0.02. 

fc =16 


hence 1 < i? < 1, 02. 

Now we have to estimate Observe that A = YlT=i — e“(2*:^-2fc(i-a;))A j 

is negative for | < a; < 1 because all terms of the series are negative. We will estimate the 
following positive quantity 


1 — a; 


J22k 


^—{ 2 k^— 2 k(l—x))/t _ ^—{ 2 k^-\- 2 k(l—x))/t 


2ke 


-2k^lt 


^ 2 k(l—x)/t _ ^— 2 k{l—x)/t 

1 — X 


1 — X 


4 2 fc^^- 2 fcVf Sinh( 2 fc(l -a;)/t) 
^ t 2k{l-x)/t 


For hxed t > 0 and k = 1,2,3,... function g{x) = decreasing for | < a; < 1 

hence its maximal value is attained for x = ^ and is equal to Thus 


1 — a; 


^ t 2k/t 




k=l V n =16 / 

Finally, using the estimates —0,15 < < 0 and 1 < il < 1, 02, we get the desired result: for 

0 < t < I and I < a; < 1 


0,85 < - A + B < 1,02, 

1 — a; 



SUPREMUM DISTRIBUTION OF BESSEL PROCESS OF DRIFTING BROWNIAN MOTION 


9 


which, for ^ < x < 1, implies the following 
0,85 < 




< 1 , 02 . 


\ sinh/. ^^2i/2_(i_3;)2/(2t) 

V > sinh(/ia;) ^ 

If we want to write factor instead of (1 — x) in the denominator, we need to multiply the 

right-hand side constant by 3/2, because for 0 < f < 1 and 1 < x < 1 we have | < ^ < 1. 
This gives the estimate for t | with constants 0,85 and 1,53. Finally, taking into account 
estimates for 0 < a: < 1 and for 1 < a: < 1 we get the desired estimate, which ends the proof 
of the theorem. □ 


Now the estimate for t>\. 
Theorem 8. For t >\ and 0 < a: < 1 


0,8 < 




< 1 , 2 . 


TT sinfTra:) 

Proof. For t | we use formula fl5.13|) and the following inequality: for 0 < a; < 1 there holds 
I sin(/c7ra;)| < fcsin(7rx). The hrst term of the series l)^’''^fcsin(A;7ra;)e“^^’^^*/^ is much 

larger than the sum of the absolute values of all the rest: 

OO CXD 


sm TTx e 


k=2 


-TT^tl2 


k=2 


E 


ne 


-(n-l)7r2/8 ^ 


e-^V4('4e^V8 _ 3) 

^ ^ sm(7rx)e < 0,2 sm(7rx)e 


n=4 


(e-V8 _ 1)2 

This means that 

OO 

0,8sin(7rx)e-"'*/2 < ^(-1)^+1^; sin(A;7ra:)e-^'^'‘/2 < i, 2 sin(7ra:)e-"'*/2^ 


k=l 


hence for f > 1/4 and 0 < x < 1 

0,87rsin(7rx)^^J^i^e-'^'*/2-,r2'/2 <7i(t,x) < 1, 27r sin(7rx) ^ , ^-r^tl 2 -.H/ 2 ^ 
smh(/ix) smh(/ix) 

which ends the proof. 

Now, because for all x G (0,1) there holds < dvr, the above inequality implies 

the following: 

l\t,x) 


□ 


(5.17) 


0, 87r^ < 


< 4, Sn. 


\ sinh^t „-^H/2-ttH/2 

Observe, that the above inequality differs from that given in Theorem [7] for the case 0 < x < 1: 
it does not contain factor in the denominator and instead of a factor 

it has a factor e~'^ If we want to have one estimate for all f > 0 and 0 < x < 1, we must 
add such factors. Multiplying the estimating function from Theorem [7] by e“^ G (e“^ 1), 

we must multiply the constant on the left-hand side by e“^ ~ 0, 29. This operation changes 

constant \ from Theorem [7] to ~ 0,0728 > 0,07. 


by 


On the other hand, the estimating function from the denominator in fl5.17p must be multiplied 

+ ^-(l-xf/(2t)_ 

\/^{t + x)t^/^ 
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But for t i and 0 < a; < 1 the above function is greater than l/\/^ 0,3989... > 0,39 
and less than ~ 4,92... < 5 hence constant 4, Svr from 05.171) must be changed to 

5 • 4, Stt = 24?! and constant 0, from fl5.17p must be changed to 0, Svr^/> 3,1499 > tt. 
In this way theorems [7] and [H] together imply the following: 


Corollary 1. For all t > 0 and 0 < a: < 1 the following inequality holds 


0,07< 


__ 

x{l-x) sinhifi) (l+t)5/2 ( 

t+x sinh(/.ix) f y 2 2t J 


< 247r < 75,4. 


It is easy to notice that by 05.14p the following scaling property holds: 




sinh(p,a;/ro) 


exp I — I 1- y 


sinh(/i) sinh(/ia;)rQ 
This, together with the above Corollary proves the following. 
Theorem 9. For /a, tq, f > 0 and 0 < x < r^ we have 

Y°{t,x) 


(5.18) 


0,07 < 


7 


t X 
2 ’ 

Tq ro 


x{ro-x) sinh(/jro) f _(rorF+]^f _ (rp-xy 

t+rpx sinh(^x) \ 2r^ 2* 


< 75,4. 


Remark. Many formulas in the book [2] are given in the language of a special function ssy{v, t) = 
^y 3 i 2 YlT=-ooi'^ — v + 2 fct)e“ 0 -^+ 2 fcb^/( 2 ?/)^ V <t (see [2], page 6fl). Observe that sst{x,ro) is 
precisely the series in our Y°it,x) and was estimated above. Using our method it is possible 
to give an estimate of sst{x,rQ) for all possible values of variables (like we did it in Theorem\E 
for Y°(t,Y) estimates for different sets oft and x, given in TheoremsY^ and\Bi are much 
more exact. For instance the proof of Theorem^ gives the following: for ro = l,0<f<| and 


0 < X < ^ there holds 


0,25 < 


SSt{x, 1) 




< 2 , 01 . 


\/27rt^ {tYx) 

Now we estimate the density of the transition probability of the killed process. 
Theorem 10. For fixed ro > 0, all x,y ^ (0, ro) and t > 0 there holds 


p'"°{t;x,y) 


jrl + tff^^ flA — 

Tq sinh(/aa;) sinh(/ai/)\/t V t 

{rofiff + Y, (x-y) 


1 A 


(ro -x)(ro -y) 


t 


X exp 


2rl 


t — 


2t 


The constants in the above estimate can be taken ci = 0,0029 and C 2 = 2413. 

Proof. Observe that 

1 \\ ^ f t X y 


(5.19) V" (i; X, y) = -P I - 


ffft 


1 - 


P 


2 ’ ’ 
r( ro ro 


sinh(/ax) sinh(/ai/)ro 

thus it is sufficient to consider only the case when ro = 1. Dehne the following function 

1 f (w + 2ky 

(w + 2A:)exp f-■ ^ 

k=—oo 


^ (i« + 2a)exp 


2t 


w G 
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Note that the function 554 ( 0 , w), mentioned in remark to Theorem [HI is given by the same 
series, only the range of its argument w is different. By the formula fl5.15p we have 

rx+y 

P^{t;x,y) =——— . —r/ \{w)dw. 

smh{nx) smh(/i?/) J^_y 

From the dehnition of X{w) we can deduce the following properties of the function: 

X{—w) = —X{w), A(1 + w) = —A(1 — w), 

which imply 

g-M^i/2 pl-ll-x-yl 


= - 


smh{yx) sinh{yy) 


X{w)dw. 


Observe that \x — y\ < 1 — \1 — x — y\ ioi x, y & (0,1). Since 7 ^(f, x) = ^™sinh'(/^x) ' ~ 

(cf. I5.14IB . we get by virtue of Corollary 1: 

(5.20) 

Consequently 


x{l - x) {1 + 

-P7^«P 


TT^ (l-a;)2 

' 2 * 2t 


p^{t;x,y) 


e-('^'+^')*/2(l +t)5/2 

sinh(/ia;) smh.{yy)t^/^ J\x-y\ t + {I — w) 


"l^^dw. 


Now, substituting w = we get 

1 — tn = 1 — \/l — vt = 

hence we have 


vt 


1 + ^/l — vt 


and dw = 


—t 


\/l — ut’ 


(5.21) 


-tv < 1 — w < tv. 


Combining this with Lemma [T] from the Appendix, we obtain 


p^{t]x,y) . — —rTTT^ 




sinh(/ia;) sinh(/ii/)fh2 

sinh(/ia;) sinh(/ii/)fh2 
1 - {x-yf 


U-{l-\l-x-y\P 1 + U 


'’/‘^dv 


X 1 A 


1 A 


(1 - \l-x-y\f - {x-yf 


We rewrite the expression in the last parentheses as follows 

x-y\f - {x-yf = 


Axy y -\- X <1, 

4(1 — a;)(l — ?/) y + x>l, 


= A[xy A {{1 - x){l - y))]. 


Moreover, it is easy to check, that 1 < < A ioi t) < x < 1 — y < 1. Since f{x,y) = 

1 — {x — y f = /(I — x,l — y), we get for all x, y G (0, 1) 

l-{x- yf 


1 < 


{xy) V ((1 -a:)(l -?/)) 


< 4. 
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This implies that for all x, y ^ (0, 1) 


A . {l-\l- x-y\f-{x-yr 


1 - 


t 


t 

^ ^ {xy) V {I - x){l - y)) ^ ^ {xy) A ((1 - a:)(l - y)) 


and similarly 


1 A 


l-{x- yf 


^ ^ (1 - \l-x-y\Y - {x-yy \ ^ 


1 A 


(1 -a:)(l -1/) 


Theorem [9]imp lies that in the estimate (15.20^ of A(x) we have constants ci = 0, 07 and C 2 = 75,4. 
Using this, the above inequalities, constants from Lemma 1 and inequality fl5.2ip . we get the 
constants in the thesis of the theorem. □ 


Recall that Mt = supo< 5 <t 2’^ and denote m{t,x,y) = P*(Mt G dy). An estimate of this 
density is the most complicated. In the proof we will use three elementary lemmas, which will 
be proved in the Appendix. 

Theorem 11. For all 0 < x < y and t > 0 the following estimate holds: 


m{t; X, y) ~ 


x{y - x) sinh(/i?/) (l + ^ + ^ + \/^/^) 

yH sinh(/ix) {y'^ + t^iimY + 1)) {t + yx) 





hy^+h(yM)^+i)) 


Proof. By the strong Markov property (SMP) we have for y' > y > x 


¥^{M,e{y,y')) 


SMP 


{Py < t, Ty' > t) 

[Py < t] E^ {Pyl > t — Py)] 
nt noo 


7 ^(m,x) / {v,y) dv du. 


• t—U 


By fl5.18|) we obtain 


P"(M,e(i/,i/')) 


x{y - x)y{y' - y) sinh(pi/') 


{y'vY 


sinh(/ia;) 


7* {y-^) 

X / —X_ i _ P ^ 5 ” 

Jq {u + yx)u^/'^ 

■=F{t,x-,y,y’). 


I - - , . „,„ e 272 dvdu 

t-u [v + y'yjv'^/^ 
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Since m{t; x, y) is continuous (it is obtained by differentiating formula fl3.9p with respect to tq), 
then m{t] x, y) = hence 


u 1 ~ r 2/') 

m{r,x,y) ~ hm- 

y'^y y' - y 

x{y — x) sinh(/i?/) 


y"^ sinh(/ia;) 


X 


(l/^ + 1 '°° _ (!£iid+ 


[U 


yx)u^/‘^ 


e 2 b 


' t—U 


y3/2 


e ^ dv du 


x{y — x) sinh(/i?/) -(ynf+ilt 


y"^ sinh(/ra;) 


e ‘^y 


X 


(l/^ + (y-^f (y"^ + t — U + _ (yyf+- 


If) {u + yx)u^/‘^ 


e 


/o {t — u + 


e 2 P dy du. 


We apply Lemma 2 from Appendix with a = y‘^,b = t — u, c = to the inner integral 


and get 


2 , _2 rt 


m{t;x,y) 


x{y — x) sinh(/ii/) -IsEil+ilt 


(y-^) 


y^ sinh(/ia;) 


e ^y 


{y"^ + u)^^‘^{y‘^ + t — u)^^‘^e 

If) s/t - u{u + yx)u‘^/‘^ {y^ + (t - u){{yiJ,y + 1)) 


du. 


We split the above-given integral into two parts = h + h and substitute w = 

(u “ f)’ ^ ~ (u “ i) -^ 1 ) -^ 2 , respectively. Before substitution we use the following 
estimates: t—u ~ t for m G (O, |) and m ~ t for m G (|, t). Additionally, the second substitution 
gives t — u = 2t^w/{2tw + {y — x)'^). Consequently 


+ + ^)) k {u + yx)u^/‘^ “ 

/ N 5/2 

^^4/^2 ^ ^)3/2g-(2,-x)2/t .00 \^ + {y-xf{\ + ^)\ e""' 

- / -- - - dw 

Vix{y - x) ( 1/2 + tiiyyy + 1)) Jo ^ ij^y (w + {y - x)2(l + ^)) 

(i/2 + t)5A A (i/2 + t_^)3/2e-i^ 

{t + yx)t^/‘^ Jt/2 \/t-u (i/2 + {t- u){{yyY + 1)) 

2(i/2 + tfo-iy-xfyt r{y-x?/2t {w + ^ 

(y-x)Ht+yxw+w+i)l 


Now we apply Lemma 3 to the integral in the estimate of A with a = (i/ — a;)^(j + ^), 
b = ^ c = {y — xY{j + and Lemma 4 to the integral in the estimate of A with 


a = 


2 t ’ 


b 




{y - x)W ^ {y - xf 1 ^ {y - xf 

2t{y‘^ + t) t 2 ( 1 -|-^) t 
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SO all assumptions of Lemma 4 are satisfied. Hence 


h + h 






+ 


Vi{y^ + t{{yy)^ + 1)) x{y - x) \{l+ (^^)(^) I + {y - xf{^) 

(i/^ + 


{y - xfit + yx){y‘^ + t{{yy)‘^ + 1)) 


X 


X I 1A + V^(i/^ + ^(W^ + i))(i/ - a;) V 


2t 




I / \ r» / 'ti'T’-l—i- \ I 


\/^ ( 2 /^ + f( + 1)) (f + 2 /a;) V Viy^{l + ^^ 7 ^) 1 + ( 2 / - a;)2(^) 


+ ( 1 A {y'^ + tf^^Vi ^ (2/^ + t)^/y^ + t{{yyy + 1) 


2f 


( 2 / - ^ /i I _ 


Let us denote the expression in the last brackets by J^{y, x, t, y). To prove the theorem we need 
to show that 


(5.22) 


J{y,x,t,y) 


(?/^ + t) 1 + ^ + ^ + Viy 
yVi 


1 + 


(y-x)‘^y^ 


t{y^+t{{yy)^+l)) 


Assume > 1. Then we have y'^ > t, which implies 


(^1^2 _|_ ^^5/2g-(i/-a:)2/2t 


< 2 - 


(l/2 + t)3/2 


< 2 


Viy^il + ^ 2 / 2 -\/ i 7 ^ 


1 + (» - A^l'gr) 


2/ 


( 2 /^ + ^)V2/^ + f((2//^)^ +^) 

( 2 /^ + f)V2/^ + ^(( 2 //^)^ +~) 

\/f2/S/^ 


{y—x^y^ 

t{y^+t({yyy^+l)) 


^ ^ {y-xY \ {y^ + tf/^^Vi ^ ^ ^ ( 2 /^ + ^)V2/^ + ^((2//a)^ + 1) 


2t 


{y - xfy^ 


Viy^^Vf 


Viy^\ I + 


(j/—a:)2j/2 

i{y'^+ii{yfJ-)'^+^)) 


SO we can estimate J'{y, x,t, y) from the above by its last component. But for x,y,z > 0 
there holds -^/x^ + i/2 -p ^2 ^ x + y + z, because I 2 and h norms are equivalent in M^. Hence 

Vy Ri 1 + ^/ty + ^ and because ( 1 / — x)^ > t, then fl5.22p follows. 

In the case < 1 we have t + xy ^ t + y'^ and also (recall that 0 < x < y) 


1 + ( 2 / - x) 


^yx + t _y , X , {y-xf 


yxt 


—-^-h 

X y 


i^y. 

X 
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Hence 


J{y X t ^ t + y^ ^ {y^ + ^ {y^ + t)^yy^ + t{{yyy + 1) 


Viy'^ 


y{y - x) 




Viy^ 


y“ + t ( 2 {y^ + ty/^ ^ y/t ^ \/y‘^ + t{{yyY + 1 ) 


y\/i 


yO 


y-x 


y 




yVi 


y. 


y-x 


y 


which again is eqnivalent to fl5.22p . 


□ 


6. Appendix 


Here we gathered fonr lemmas which were nsed in the estimates carried ont in the previons 
section. 


Lemma 1. For 0 < a < b we have 


fb w^w/2^^ 
Ja l-\-w 


1 

12 - iV2(iA6)(l A(6-a)) 


< 2. 


Proof. Let us denote X(a, h) = j^e'^^'^dw and let a,b < 1. For w G (0, 1) the function 

fi'^) = is decreasing hence for such w there holds ^ < f{w) < 1. Using these inequalities 

we get 




w dw < 


w 


1 + tc 


e'^^‘^dw< / wdw = b‘^ — of = {b —a){b +a). 


Since 0<a<6<l, we have b <b + a <2b and this implies 


— a) < -^b{b — a) < X(a, b) < f wdw < 2b{b — a) < 2e^^‘^b{b — a). 

^ ^ J a 

Now let 6 > 1. We use inequality 1 — e~^ < 1 A x, x > 0, and get 

X(a, b) < f e^/^dw = (l - < e^/^ (1 A (6 - a)). 

J a 

For X > 0 it holds that 1 —e“* > > ^(1 Ax). Using this, the fact that the function x —)■ 


is increasing and (a + 6)/2 > 1/2, we obtain 


X(a, b) > 


w 


I g + b \ W 


-e^'^dw > 


a+6 


1 + 4 ^ Ja+b 


e->/^dw > Y^2e'/2 (1 - 


>_eU2(iA(6-a)), 
which ends the proof. 

Lemma 2. For a,b,c> 0 and ac> 1 we have 

(6,23) + 6 + 


□ 


’ dv 


_|_ 5)3/2 


(6 + n)3/2 " ~ y/b{l +be) 

Proof. Let us denote the integral in the thesis by /(a, b, c). Substituting v = ^ we obtain 

/•OO 


( 6 . 24 ) 


T( -L \ "I ((fl + b)c + s)^/^ _ 
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For be > 1 we get 

/(a, b, c) = 


1 /a + 6(1 + s/((a + 6)c))^/^ _ 1 /a + b 


- f 9 5^ — 

(1 + s/(6c))3/2 C 


3/2 


which is equivalent to 16.231 For 6c < 1 we split the integral in fl6.24p as follows: “ /q^ +“ /i°° = 
Ji + / 2 . Then 

I 2 ~ Vc(a + 6)^^^. 

In Ji we substitute s = bet and get 


hence 


h 


h + h 


-l/6c ^ 3/2 


(1 + t)3/2 


: dt K. b 






3 / 2 ^ 


From the assumptions be <1 and ac > 1 it follows that 


be 


a + b 


3/2 


> (be 


fn ( 


a + b 


\ b 


3/2 


> 1, 


which ends the proof of the lemma. 

Lemma 3. Let a,e > b > 0 and a < 2 if b < 1. Then we have 


□ 


+ a)5/2 


1 

+ 


(6.25) / . ^ - re-^dw 

Jq {w + by{w + e) b{b + l)e 1 + c 

Proof. Denote the integral in the thesis by F{a,b,e). It is clear that for 6 > 1 we have 
F(a, b, e) ~ which is equivalent to fl6.25p . For 6 < 1 we split the integral into two parts 


F{a, b, e) ^ 


If) {w + by{w + e) 


e ^dw + 


w 


, 5/2 


If) {w + by{w + e) 


■e-^dw. 


The last integral can be approximated by which is shown below: 




W 


1 + c Jo {w + ly w + e 


e-'^dw < 


w 


, 5/2 


< 


If) {w + b)‘^{w + e) 

POO 


e ^dw 


w 


0 ^/ww + c 


-e ^dw 


1 + c 


Moreover, 


Iq {w + b)‘^{w + e) 


.e-^dw = j +j 

_ 1 
be Jo 
1 1 


3 ^'^du 


■ + 




(1 + 11 ) 2(1 +-m) Jf) {w+ b + iy{w + e+1) 


be 1 + c ’ 

Hence for 6 < 1 (in this case a < 2 by assumption) we obtain 

1 n5/2 


F(a, b, c) + 


6c~^l + cy~*”l + c be ~^l + c’ 


1 

^-h 


which is again equivalent to fl6.25p . 


□ 





































SUPREMUM DISTRIBUTION OF BESSEL PROCESS OF DRIFTING BROWNIAN MOTION 


17 


Lemma 4. Let a > b > c > 0. Assume also that (a) a ~ 6 z/a > 1. Then we have 
(6.26) 


iw + '^dw , . 

^ ’ ^(lAa) + 


0 + c) Vc(l + c) ‘ 

Proof. Let us denote the integral by I {a, b, c). Assume now that c > 1. Then we have a,b > 1 
and get 

63/2 ra 02 


I {a, 6, c) = 


+ 1 ) 


It is equivalent to the thesis because > 6^/^ > 1. In the case c < 1 we substitute 

w = cu: 


(6.27) 


/(a, b,c) = c 


ra/c (^, I ^—cu 




-du. 


Iq y/u{u + 1) 

For a > 1 we use the assumption (a). It follows that u + - ~ - for u G [0, -]. Hence 


I {a, b, c) 


b\3/2 ^a/c g-^ 


- 


-du 


02 


sj Jo y/uiu+iy^ y/c' 

which, as before, is equivalent to the thesis. For a < 1 the formula fl6.27p gives us 

udu 63/2 ra/c q-uc 0 


I {a, 6, c) 


'0 

The last integral is bounded by from below and by ^J+i) fr°^ above. Moreover 

/■“fr udu 


n + 1 


+ 


\/u{u + 1) 


du. 


/a / a 

= c(- ln(l + - 

M + 1 Vc V c 


which ends the proof. 


□ 
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